Fusion hindrance at deep sub-barrier energies 
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Fusion of two heavy nuclei may be initiated at larger distances (before contact point is reached) 
due to neck formation or due to strong neutron flow between the two colliding nuclei leading to 
an elongated di-nuclear configuration. When injected into the adiabatic fusion trajectory, these 
configurations undergo further up hill diffusion to achieve fusion, thus causing fusion hindrance 
at deep sub-barrier energies. Coupled channel calculations based on this new concept is able to 
explain the measured steep falloff in fusion cross sections and associated anomalies at deep sub- 
barrier energies. Further, it is also shown that fusion initiated by neutron flow followed by diffusion 
can provide a natural explanation to fusion of radioactive ions charcterirzed by extended neutron 
halos, a topic which is highly controversial. 

PACS numbers: PACS numbers:25.70.Jj24.10.Eq,25.60.-t,25.70.Gh 



I. INTRODUCTION 



The study of near and sub-barrier fusion of heavy and 
exotic ions have generated a considerable amount of in- 
terest and anomalies over the last several years. In the 
simplest picture, fusion reaction between two ions at sub- 
barrier energies is governed by the penetration through 
the Columb barrier followed by an absorption inside the 
barrier. Early anomalies refer to the failure of this sim- 
ple barrier penetration model (BPM) which could not ex- 
plain the large enhancement found in the measurement 
of fusion cross sections for all heavy ion systems particu- 
larly at sub-barrier energies 0, Q. Since then, several 
theoretical models have been proposed to account for 
this large sub-barrier enhancement. The most success- 
ful model that has emerged out of these studies over the 
past several years is the description of fusion within a 
coupled channel framework where the presence of cou- 
plings to various low lying inelastic and transfer chan- 
nels are treated explicitly. Under certain approximation, 
the channel couplings can be interpreted as a phenom- 
ena by which the colliding nuclei encounter a distribu- 
tion of barriers, one or more of which having height less 
than the original Coulomb barrier, thus giving large fu- 
sion enhancement. Therefore, a better understanding of 
the fusion process followed through the distribution of 
barriers which can be derived experimentally from the 
second derivative of the product of the fusion cross sec- 
tion with energy Q. In most experiments, fusion cross 
sections have been measured down to the 1 — O.lmb level 
and the coupled channel calculations have been quite suc- 
cessful in reproducing the general trends of the measured 
yields. While efforts are still being put to improve upon 
the techniques of these calculations to understand the 
nature of fusion enhancement and the associated bar- 
rier distribution, some of the recent measurements at 
extreme sub-barrier energies [U, [j| (cross section below 
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the lib level for systems like 64 Ni + 64 Ni, 58 Ni + 58 Ni, 
60 Ni + 89 Y) have brought out new surprises which show 
supression in fusion cross section with respect to the same 
coupled channel calculations that explains the enhance- 
ment at sub-barrier energies. The most striking feature 
of these measurements are the logarithimic derivative 
L(E) = d(lntjE) / dE which continues to increase steeply 
with decreasing energies where as all the realistic model 
calculations including that of coupled channel results pre- 
dict a slow rise in L(E)oi eventually tends to saturate to 
a constant value @, 0, Q . All these models also fail to ex- 
plain the behavior of 5-factor which shows the presence 
of a maximum indicating a strong hindrance to fusion at 
deep sub-barrier energies [H, 0] ■ 

Recently, these anomalies have been explained success- 
fully by using a drastically modified nuclear potential in 
the barrier region that results in a shawlloer potential 
pocket as wll as a thicker potential barrier 0, [l(| • While 
the thicker barrier is responsible for the steep falloff of 
the fusion cross secion at deep sub-barrier energies, the 
shallower potential pocket prevents higher partial waves 
to fuse. Although the proposed concept is highly illumi- 
nating and has several new features [lfj that requires fur- 
ther experimental investigations, the model suffers from 
one technical draw back due to the existence of a higher 
cut-off energy below which the system will not be able 
to fuse. For example, in case of 64 Ni + 64 Ni system, 
the minimum value of the potenial pocket V m i n rs 85 
MeV which is about 9 MeV below the Coulomb bar- 
rier. This would imply that fusion cross section Of will 
vanish for E cm < V m i n although compound nucleus for- 
mation is still possible as long as (E cm + Q) is positive 
i.e. compound nucleus can in principle be formed down 
to a threshold energy of E cm « 49 MeV. The experi- 
mental measurements also do not suggest a sharp cut-off 
in <jf for E cm < V m i n . Therefore, without putting any 
stringent requirement to modify the nuclear potential at 
the interior of the Coulomb barrier, we propose an alter- 
nate model, yet very sinple to implement within a cou- 
pled channel frame work, by allowing fusion to initiate 
through neutron flow or through neck formation at a dis- 
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FIG. 1: (Upper panel) The potential V(R) versus inter-nuclar 
distance R for 58 Ni + 58 Ni and ie O + 208 Pb systems. Lower 
panel shows the two centre shell model potentials for indi- 
vidual nuclei (dashed-dot curves in red) and merged neutron 
potentials (solid curves) for 58 iVi+ 58 Ni and 1(i O+ 208 Pb sys- 
tems at distances of 12.5 fm and 14.1 fm respectively. In the 
lower panel, x-axis has been shifted by 50 fm for clarity for 

16 Q + 208 p h system 



tance much before the Coulomb barier is reached. Since 
these di-nuclear complexes (two nuclei connected by a 
thin neck) are formed at large distances, further up hill 
diffusion over an adiabatic fusion trajectory is required 
in order to achieve fusion. This new concept not only 
able to explain recent experimental anomalies at deep 
sub-barrier energies, but also it accounts for the general 
trends of fusion process for several heavy ion systems 
over a wide range of energies starting from above to deep 
sub-barrier energeis. Contrary to the general expectation 
that the two nuclei will fuse by tunneling through the 
Coulomb barrier, the fusion in the present model is pro- 
posed to be initiated by neutron flow at large distances 
followed by disffusion over the Coulomb barrier. Inter- 
estingly, apart from heavy ion systems, this model also 
explains fusion phenomena of radio active nuclei (having 
extended neutron halos) in a natural way, a topic which 
is highly controversial pH El El EI EE EH . 

It is the normal practice to assume that the collid- 
ing nuclei will fuse tunneling through the barrier when 
E < Vb where Vf, is the height of the Coulomb barrier. 
Fig. Q] (upper panel) shows the plots of total potential 



as a function of inter nuclear distance R for two typical 
( 5S Ni + 58 Ni and le O + 208 Pb) systems. It may be men- 
tioned here that when the colliding nuclei approach each 
other, this Coulomb barrier provides a respulsive force 
to the proton clouds only and does not prevent neutrons 
to flow from one nuclei to other. In fact, it was argued 
by Stelson [17| earlier that fusion may be initiated at 
large distances when it is possible for neutrons to flow 
between the colliding nuclei. This we investigate for a 
few systems by using the neutron shell model potentials 
as was used in [l7|. We assume a Woods-Saxon shape 
with V = Vol {I + exp[(R - Ro)/a}}, where V = -50 
MeV, Ro = 1.24A 1 / 3 fm, and a = 0.68 fm. Fig. Q] (lower 
panel) shows the example of the neutron potentials for 
64 Ni + 64 Ni and w O + 208 Pb systems centered on each 
nuclei (dashed-dot curves in red) as well as the merged 
neutron potentials (solid curves) at a distances of 12.5 
fm and 14.1 fm respectively. We take one half of the two 
neutron separation energy S(2n)/2 as an approximate 
measure of the binding energy of the valence neutrons. 
We estimate 5(2n)/2 w 11.2 MeV and S(2n)/2 w 7.0 
MeV for above two systems. In the figure, the inter nu- 
clear separation has been adjusted to 12.5 fm and 14.1 
fm respectively so that the barrier of the merged poten- 
tial is deeper by the corersponding 5(2n/2) values. This 
is the situation when there is no barrier and the neu- 
trons can flow freely. This distance R t is the thrshold 
distance (marked by the arrow in the upper figure) be- 
low which there is no barrier for neutrons to flow. In 
terms of energy, there is a threshold energy E = Vt (the 
total potential at R = R t ) above which neutron flows 
freely and for E < Vt, neutrons need to tunnel through. 
In [J] , from the experimental measureemnts, the authors 
have derived a thershold energy E t for a number of sys- 
tems. The E t is the threshold enegy below which the 
fusion cross sections falls off steeply (S factor shows a 
maxima). Fig. [2] shows the plot of E t (which are derived 
from the experimenal data [J] and listed in second column 
of table I) and Vt for the same system estimated as de- 
scribed above. Interestingly, a strong correlation is found 
between the two threshold energies E t V t (V t ~ 0.95E t ), 
which is an indication that neutrons flow may be acting 
as a door- way state for fusion [It], EH • Macroscopically, 
it is equivalent to say that neutron flow provides a force 
strong enough to overcome the Coulomb replusion thus 
promoting neck formation and forming a di-nuclear com- 
plex at a distance away from the barrier position R^. 
This neutron flow is rather sequencial and distinctly dif- 
ferent from one or two neutron transfer which is generally 
treated as simultaneous coupling in the conventional cou- 
pled channel calculation [lj| . These di-nuclear complexes 
can be best described as a swarm of Brownian particles 
suspended in a fluid at temperature T instead of consid- 
ering them as two separate entities approaching towords 
each other at a relative centre of mass energy E cm . Al- 
though, some of the swarm will be sliding down due to 
replusive Coulomb force, because the swarm's width in- 
creases with time, the particles will diffuse up-hill and a 
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FIG. 2: Threshold potential Vt (sixth column of thable I) 
versus threshold energy Et (fifth column of table I). The solid 
line is a linear fit to the data points (Vt w 0.95-Et). 



fraction will be able to overcome the barrier and achieve 
fusion. As will shown later, fusion achieved through this 
diffusion process shows many interesting features that 
have been observed at deep sub-barrier energies. 

The paper is organized as follows. After this brief in- 
troduction, in section II, we derive an expression for the 
one diemsnional diffusion model of fusion and compare 
the results with the predictions of the standard barrier 
penetration model. In section III, we propose a cou- 
pled channel formalism where fusion is estimated as a 
weighted sum over the distribution of barriers, but with 
a transmission factor estimated using the diffusion pro- 
cess. The results of the coupled channel calculations are 
compared with the experimental measurements. In sec- 
tion IV, we extend this formalism to explain experimental 
fusion cross sections of radio-active halo nuclei. Finally, 
the conclusions are presented in section V. 

II. ONE DIMENSIONAL FUSION MODEL 

A. A diffusion model for fusion 

First, we consider an one dimensional diffusion model 
under the influence of a linear driving force (F = 2bR) 
corresponding to a parabolic potential, 

V ad = B-b{R-Rff , (1) 



TABLE I: The threshold energy and other relevant parame- 
ters for various systems. The strength of the nuclear potential 
has been adjusted to reproduced the same Columb barrier as 
used in Ref. [4J and are listed in the second column. The corre- 
sponding R), values are listed in the fourth column. The fifth 
column lists the threshold energy E t at which the <x/ falls off 
steeply. The last two columns show the Vt and Rt values de- 
rived using the method as described in the text. Note that 
for all the systems, Rt > Rt as expected. 



System 


v b 


Rb 


S'2n/2 


Et 


V 


Rt 


b *Ni + bs Ni 


101.2 


10.3 


11.2 


93.9 


90.0 


11.3 


w Ni + «y Y 


132.5 


11.1 


10.2, 10.4 


123.o 


117.1 


13.4 




187.7 


11.2 


10.6, 10.4 


170.8 


159.0 


14.1 


iM Zr Zr 


189.6 


11.3 


10.6, 7.9 


170.8 


157.7 


14.6 


m Zr + yu Zr 


192.5 


11.2 


10.6 


175.2 


163.1 


14.1 


Ti + Pb 


203.6 


12.1 


9.5, 7.1 


181.6 


167.5 


15.5 


M Ni + M Ni 


98.0 


10.7 


8.2 


89.2 


84.7 


13.3 


i« + ^us pb 


76.3 


11.7 


14.4, 7.1 


70.1 


66.9 


14.1 


m Ca + yu Zr 


94.3 


10.4 


14.5, 7.9 


94.3 


87.1 


13.2 


ib O+ 144 Sm 


56.6 


10.9 


14.4, 9.6 


56.6 


55.6 


12.8 


w F + '2m pb 


86.0 


11.6 


9.8, 7.1 


75.1 


74.7 


14.2 


10 O + 10 " Sm 


60.6 


11.0 


14.4, 6.9 


52.4 


53.2 


13.4 


Wp + 2UH pb 


86.0 


11.6 


9.8, 7.1 


75.1 


74.7 


14.2 


m Ar + 144 Sm 


133.5 


11.1 


8.2, 9.6 


116.6 


113.8 


14.1 


4,1 Ar + lb4 Sm 


131.4 


11.3 


8.2, 6.9 


108.8 


110.7 


14.5 


4il Ar+ n ' 2 Sn 


110.6 


10.9 


8.2, 8.5 


96.6 


95.2 


13.6 


« b Kr Mo 


179.9 


11.1 


8.5, 11.4 


160.4 


151.0 


14.4 


* b Kr + 1W Ru 


185.9 


11.5 


8.5, 7.6 


162.4 


153.9 


13.4 


w s+ iw pd 


90.3 


11.0 


8.4, 7.5 


79.7 


77.8 


13.5 



where B is the barrier height at R = R% . Note that 
we have introduced the notation ad to indicate the adi- 
abatic nature of the potential and to make it distinguish 
from the sudden (frozen density approximation) potential 
which is generally used to describe the fusion trajectory. 
For example, the potenial shown in fig. [T] is commonly 
refered as sudden potential and is more appropriate to 
describe reaction dynamics when the relative motion be- 
tween two colliding nuclei is fast enough as compared to 
other intrinsic degrees of freedom. Under this sudden 
or frozen density approximation, the potential depends 
only on the relative separation between the two collid- 
ing nuclei. However, due to neutron flow, if the colliding 
system forms a di- nuclear complex at R > R^ d , the po- 
tential is no longer a function of relative separation alone, 
but also it depends on the neck degree of freedom (two 
nuclei connected by a thin neck). The resulting fusion 
trajectory will be more adiabatic in nature and needs to 
be obtained by optimizing a set of nuclear shape param- 
eters, a procedure what is generally adopted to obtain 
a conditional fission trajectory. However, we expect this 
adiabatic fusion trajectory to be quite different from a 
true fission trajectory at least in the exterior region as in 
case of fusion, the neck growth is not as rigorous as in the 
case of fission process. Although it will be interesting to 
estimate the fusion trajectory by parametrizing the po- 
tential in terms of a set of shape parameters that can 
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generate various di- nuclear configurations, in the present 
work, we keep the formulation simple by consideaing a 
parabolic form of potential as given by Eq.{T]). Thus, we 
denote Vb for sudden barrier height where as B represents 
the height of the adiabatic potential barrier. 

Defining W(R, t) as the probability of finding a Brow- 
nian particle at position R and time t, the equation de- 
scribing the drift and spreading of W(R,t) is given by 
Smoluchowski diffusion equation which is a special case 
of a Fokker-Plank equation [2(| HH , 



d t W(R,t) = Dd R {d R - F/T)W(R,t), 



(2) 



where 9r is the partial derivative with respect to R, T 
is temperature, F — 2bR is the driving force and D is the 
diffusion coefficient. The distribution W(R, t) satisfies 
the initial condition at R = Ro and t = to , 



W(R,t )=S(R-Ro). 



(3) 



The above equation has a stnadard solution [20j, |2l| given 
by, 



W(x,t) 



1 



exp 



(x- < x >y 

2e 



where x = R - Rf, 



< x >— Ro exp 



2bD 



(t - t ) 



and 



AbD . 

exp \ -r=-(i - t ) I - 1 



(4) 



(5) 



(6) 



Note that both width £ and average position < x > in- 
crease with time as expected, but the ratio < x > /\/2~£ 
approaches a constant value of Ro\Jb/T as t — > oo. 
The fusion probability can be estimated by integrating 
W(x, t) from x = — oo to x = (i.e. R — — oo to R = R},) 
and taking the limit t — * oo, which turns out to be, 



T d = \[\Terf{Z)], 
where erf is the error function and the argument 



(7) 
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FIG. 3: The upper panel shows actual sudden potenial V(R) 
as a function of R for 58 Ni + 58 Ni system (solid curve). The 
dashed-dot curve (in red) is obtained with a fit parameter 
n = 1.3. The dotted curve (in blue) represents a parabolic 
potential for n = 2.0. The parameter c = 4.0 is chosen to get 
the fit as close as to the actual potential. The lower panel 
shows the similar plot for 6 He + 209 Bi system with n = 1.1 
and c = 1.0. 



comes to the rest, developes a neck forming a di-nuclear 
complex and gets captured into a trajectory which is adi- 
abatic in nature. Depending on the injection point, the 
system will undergo further diffusion. Fusion is achieved 
by diffusing up hill of the adibatic trajectory which has 
a barrier height AB at the injection point Ro- Since we 
assume Ro as same as the classical turning point, it can 
be obtained numerically by equating E with the sudden 
potential V sud (R) at R — Ro. However, for mathemat- 
ical convenience, we parametrize this potential by using 
a more generic form, 



Z = 



! b(R -R b y 



AB 



(8) 



The minus (plus) sign corresponds to the case when the 
injection point Ro > Rb (Ro < Rb)- The AB is the 
height of the barrier seen at the injection point Rq. 

We describe the fusion of two colliding nuclei in the 
following way. Initially, the relative motion of two nu- 
clei are fast enough, the role of neutron flow or the neck 
opening may not play a significant role and the colliding 
system follows a fusion trajectory which is sudden in na- 
ture. However, at the classical turning point, the system 



V sud =V b -c(R- 



r>siu 



y 



(9) 



where i?^ ud is the position of the sudden barrier Vb- 
The fit parameter n is adjusted to reproduce the sudden 
potential in the outer region (R > i?^" d ). Using the 
above parametrization, Ro can be written as, 



Ro — R, 



sud 



Vb-E 



l/n 



(10) 



Fig. [3] shows the plots of a realistic sudden potentials 
(solid curves) for two typical 58 Ni+ 58 Ni and e He+ 209 Bi 
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systems. The dashed-dot curves (in red) are obtained 
with the fit parameters n — 1.3 and n = 1.1 respec- 
tively. It may be pointed out here that for n = 2, the 
potential is parabolic which is most commonly used in 
the literature. However, we have seen a parabolic ap- 
proximation grossly underestimate the potential at large 
distances (R > R b ud ). It is noticed that for most of the 
systems in the region R > Rl ud , the sudden potential 
can be approximated with the parameter n which lies in 
the range 1 < n < 1.3 i.e the tail of the potential more 
closely follows a linear decrease, rather than parabolic. 
It is also assumed that the injection point is same as the 
classical turning point Rq at which the relative velocity 
becomes zero. This assumption is an over simplification 
as the capture into the adiabatic trajectory due to neck 
formation may take place even before the system comes 
to complete halt i.e. even before Rq is reached. There- 
fore, n can be kept as a fit parameter. However, it is 
observed that experimental data can be reproduced well 
when n is close to unity (n < 1.0). Thus with the choice 
of n = 1.0 and using Eq. (fT0|) . the argument Z in Eq.f8|) 
can finally be written as, 



(V b - E + 5) 

dVT 



(ii) 



where d = c/Vb and 5 = c(R s b ud - Rf). The term 5 
can be droped out under the assumption R b 



Even if Ri ud ^ R b d , 5 can be taken as small energy in- 



b f ±L b 

dependent positive correction which can be absorbed by 
redifining the sudden potential V b = Vb + S. However, 
in the present work, we will ignore this S term i.e. we 
will assume that both sudden and adiabatic trajectories 
have the same barrier position Rb although the corre- 
sponding barrier heights Vb and B are different. In any 
case (with or without <5), it is interesting to note that 
although Eq.((7j) refers to the diffusion probability over 
an adibatic barrier of height B, the argument Z contains 
only the sudden barrier height Vb- This is because, in the 
present formalism, the adiabatic barrier height V AB as 
seen from the injection point Rq is proportional to the 
height (Vb — E)/d. This is an important consequence of 
our model as the diffusion probability is still governed 
by the sudden barrier Vb, adiabatic effect enters obly 
through the parameter d. Substituting y = —Z and using 
the property erf(-y) = — erf(y), Eq.((7j) can be written 
as a single equation, 



TABLE II: The various parameters Vb, Rb, huj and d used to 
fit the data for five systems considered in this work. The last 
column shows the r/o values used in S(E) for normalization. 
Note that the values of Vb and Rb used to fit the data and 
the values used in table I are different due to different nu- 
clear normalization factors. In table I, 14 and Rb are same 
as used in Ref.[|[, where as, to fit the experimantal data us- 
ing CCBPM and CCDIFFUSION, a different set of V b and 
Rb are used. For all the systems except 64 Ni + 64 Ni, we use 
same V b both for CCBPM and XXDIFFUSION calculations. 
Hoeever, for 64 Ni + 64 Ni system, we found a slightly lower 
barrier explains the data well. So we have used Vb = 94.9 
MeV for CCDIFFUSION and V b = 94.5 MeV for CCBPM. 



System 


14 (MeV) 


% (fm) 


hto(MeV) 


d 


Vo 


b *Ni + bs Ni 


99.5 


10.6 


3.9 


1.2 


69.99 


64 Ni + 64 Ni 


94.9(94.5) 


11.2 


3.7 


1.5 


75.23 


16 O + 144 Sm 


60.9 


11.0 


4.6 


1.5 


41.0 


a He + 209 Bi 


20.7 


10.8 


4.2 


5.0 


17.2 


He + 238 u 


22.8 


10.9 


4.2 


3.8 


18.5 



B. Fusion through tunneling 

We can also estimate the transmission probability 
through barrier penetration using the WKB approxima- 
tion, 



TF(E) = [l + exp(A l )}, 
where Ai is the classical action given by, 
2n 



-1/ = J r {V l {R)-E)dR, 



and 



V l (R) = V c (R) + V N (R) 



l{l + i)h 2 

2/iR 2 



(13) 



(14) 



(15) 



Under the parabolic approximation, Eq. p3l) can also be 
estimated using Hill- Wheeler expression [22j |. 



<-<>A 2 £M-e) 



-i -l 



(16) 



T d (E) = - [1 + erf(y)} , 



(12) 



Finally, we can estimate the fusion cross section from 
the relation, 



where k is the relative wave number and Ti(E) is the 
transmission probability for which either Eq. (|12p or 
Eq. (fl~3"|) can be used. Even Eq. (fl6| can also be used to 
estimate fusion cross section under parabolic approxima- 
tion. 



where y = (E — Vb) / 'd. The above equation gives only 
s— wave diffusion probablity. The I— dependent diffusion 
probability can be calculated by replacing Vb with V b l — 
V b + l(l + l)H 2 /(2fiR 2 b ). 



C. Fusion cross section 

In the following, we estimate fusion cross sections us- 
ing one dimensional barrier penetration and diffusion 
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FIG. 4: (Upper panel) Fusion cross section as a function of 
E cm for 58 Ni+ 5S Ni system calculated using one dimensional 
barrier penetration (both WKB and BPM) and diffusion mod- 
els. The barrier parameters are Vb = 99.5 MeV, Ro = 10.6 
fm and hio = 3.9 MeV for BPM calculations where as d = 1.4 
is used in diffusion calculation in order to have similar cross 
section at above barrier energies. Lower panel shows corre- 
sponding L(E) versus E cm for the same system. 



model. Under the barrier penetration model, we will 
use both WKB (Eq. (flU)) ) as well as parabolic approxi- 
mation (Eq. (|16[) ). The later approximation in the text 
is refered as BPM model which has three parameters, 
the barrier height Vb, the barrier width Hlo and barrier 
position R},. The diffusion model has also three param- 
eters, Vb, Rb and d. While Vb and Rb are taken same as 
that of BPM parameters, the d parameter is adjusted to 
reproduce the experimental measurements. The temper- 
ature r is obtained from the relation E* = aT 2 where 
E* is the excitation energy and a is the level density pa- 
rameter which is taken as (A\ + A2)/10. FigH] shows 
the fusion cross section for a typical 5S Ni + 58 Ni system 
using WKB, BPM and diffusion models. The potential 
parameters are given in table II. The dotted curve (in 
blue) is obtained using the BPM (parabolic barrier pen- 
etration model) where as the dashed-dot (in red) curve is 
an exact WKB calculation. The solid curve is obtained 
using diffusion model with the same barrier parameters 
as that of BPM but d is adjusted so as to have same fu- 
sion cross section at above barrier energies. The lower 
pannel shows the corresponding logarithimic derivatives 



FIG. 5: The S factor versus E cm (upper panel) and < 
L > versus E cm (lower panel) for 58 iVi + 58 Ni system using 
BPM, WKB and diffusion models. The potential and diffusion 
parameters are same as that of figure [4] 



which is defined as 



L(E) 



d[ln(Ea)] 
dE ' 



(18) 



It is interesting to note that L(E) in BPM model sat- 
urates to a constant value, where as WKB calculation 
shows a slow rise with decreasing energy. Contrary to 
both WKB and BPM calculations, the diffusion model 
result shows a steep rise in L(E) as energy decreases. 
This is also evident from the upper panel as fusion cross 
section in diffusion model falls off more steeply at deep- 
subbarrier energies, although both BPM and diffusion 
models give nearly same cross section at near and above 
barrier energies. Another pragmatic way to represent the 
fusion process at sub-barrier energy is the S factor which 
is defined as 

S(E) = Eae 2 ^-^, (19) 

where r\ — Z\Z^ j (hv) is the Sommerfeld parameter and 
v is the beam velocity. Fig. [5] shows the plot of S factors 
obtained using BPM, WKB and Diffusion models. For 
sack of completeness, the lower panel also shows another 
measurable quantity, the average angular momentum < 
I > as a function of energy E cm . The average < I > is 
defined as, 



< I >- 



(20) 
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FIG. 6: Fusion cross section (upper panel) and L(E) versus 
E cm (lower panel) for 58 Ni + 58 Ni and 64 iVi + 64 Ni systems. 
The potential parameters are listed in table II. The lower 
panel shows the corresponding L(E) as a function of E cm for 
the same systems. The data points are taken from 0, 0|. For 
64 iVi + 64 Ni, the X-axis is shifted by +20 MeV for clarity. 



It is interesting to note that the 5* factor in diffusion 
model shows a well defined maxima at deep sub-barrier 
energy (see Fig. [5]), a feature which is supported by ex- 
perimental measurements. On the other hand, S factor 
obtained using barrier penetration models (either WKB 
or BPM) keeps on increasing with decreasing energy. 
Similarly, the average < I > value shows saturation at 
deep sub-barrier energy where as < I > obtained using 
diffusion model decreases with energy without having any 
saturation point. As the above models are one dimen- 
sional in nature, we do not compare the results with the 
experimental measurements. We only show it to demon- 
strate that the diffusion model has certain basic features 
at deep sub-barrier energies which are consistent with 
experimental observations and can not be reproduced by 
any type of barrier penetration models with the same 
type of potential as used in the diffusion model. 



III. COUPLED CHANNEL CALCLATION 

So far, we have considered transmission either through 
diffusion or through barrier penetration only from a sin- 
gle barrier. It is now well known that channel coupling 
plays an important role in heavy ion fusion enhancement 



particularly at sub-barrier energies. Although, there are 
several ways to carry out coupled channel calculations, 
here we adopt a very simple approach as originally sug- 
gested in |23| followed by the computer code CCFUS [24| . 
In this approach, instead of a single barrier V b , the fus- 
ing system encounters a distribution of eigen barriers (or 
eigen channel a) with probability Lu a . The total fusion 
cross section is written as the weighted sum of the fusion 
cross sections over all the eigen channels a and is given 
by, 



(21) 



where uj a is the weight factor for channel a and 
tjf(V b a ,E) is the corresponding fusion cross section in 
channel a which has a barrier height V b a , In our formal- 
ism, since the system encounters both sudden and adia- 
batic trajectories, we would like to argue here that the 
inelastic coupling (vibrational or rotational) will affect 
only the sudden trajectory while the adiabatic trajectory 
by and large will remain unaffected. As an example, in 
case of static deformation of the nuclei, the fusion bar- 
rier depends on the projectle-targct orientation and the 
total fusion cross section is written as the weighted sum 
over all the angles. Microscopically, it is equivalent to 
the rotational couplings under the approximation of zero 
intrinsic excitation energies. However, if the projectile 
and target form a di-nuclear complex due to neck forma- 
tion, the resulting adiabatic trajectory will be insensitive 
to the initial deformation effect as it is obtained by op- 
timizing a set of nuclear shape parameters. Therefore, 
when the sudden tajectory is replaced by a distribution 
of eigen trajectories due to channel couplings, it has no 
direct effect on the adiabatic trajectory. However, the 
injection point which is the same as the classical turning 
point Rq will now be channel dependent and needs to be 
replaced by a distribution of injection points depending 
on the channel a. Since Ro depends on R b + (Vb — E), 
V AB a should also be replaced by (V b a — E) for each 
channel a with the same weight factor uj a . The net ef- 
fect is that the diffusion probabilty will also be different 
depending on the eigen channel a. 

In the computer code CCFUS, the transmission 
propobaility from a single barrier is used using the Hill- 
Wheeler barrier penetration factor T b given by Eg . (fTT))) . 
We can use the same CCFUS code by replacing T b by the 
diffusion probability T d as given by Eq. (fT2"|) . In fact, we 
will use both the factors T b (to be refered as CCBPM) 
and T d (to be refered as CCDIFFUSION). We consider 
2+ , 3" and 4+ inelastic states of 58 Ni and 64 Ni with 
the same coupling parameters as given in [23|. For the 
CCBPM calculations, we have as usual three parame- 
ters Vb, R b and fuv and the values used in this work are 
listed in table II. In CCDIFFUSION, we use the same 
value of Vb and Rb as CCBPM, but adjust the variable 
d to reproduce the experimental data. For 58 Ni + 58 Ni 
system, the Vb, hw , Rb and d values are 99.5 MeV, 3.9 
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MeV, 10.6 fm and 1.2 respectively. The same set of pa- 
rameters are used both in CCBPM and CCDIFFUSION 
calculations. However, for 64 Ni + 64 Ni system, we use 
V b = 94.9 MeV, H = 3.7 MeV, R b = 11.2 fm and d = 1.5 
for CCDIFFUSION calculations. For CCBPM, we have 
noticed a slight lower value of Vb gives better result. Ac- 
cordingly, we use Vb = 94.5 MeV for CCBPM where as 
all other parameters are kept same as that of CCDIFFU- 
SION. Fig.e © (upper panel) shows fusion cross sections 
for the above two systems using coupled channel calcu- 
lation with both CCBPM (dashed dot curve in red) and 
CCDIFFUSION(solid curve). The one dimensional BPM 
model without any coupling effect has also been shown 
for comparison (dotted curve in blue). As expected, the 
BPM model underestimates the fusion cross section at 
sub-barrier energies. The CCBPM accounts for the fu- 
sion enhancement at sub-barrier eneries, it over predicts 
the fusion cross section at deep sub-barrier energies. The 
important observation is that the L(E) obtained from 
the CCBPM model (see the lower panel) tends to satu- 
rate to a constant value at deep sub-barrier energies like 
the BPM model as shown in Fig. 0J However, the predic- 
tion of CCDIFFUSION explains the data quite well both 
at sub-narrier as well as at deep sub-barrier energeis. It 
also expalins the general trend of L(E) which keeps on 
increasing with decreasing energy. 

Fig. [7] shows the S(E) factors (upper panel) for the 
above two systems. The normalisation factors 770 are 
69.99 for 58 Ni + 58 Ni system and 75.23 for 64 Ni + 64 Ni 
system. As shown in the figure, the predictions of CCD- 
IFFUSION are in excellent agreement with the observed 
S factors where as CCBPM keeps on increasing without 
showing the maxima. The lower panel shows the < I > as 
a function of energy. In addition to 58 Ni + 58 Ni system, 
we have also estimated < I > distribution for le O+ 144 Sm 
system. In the coupling scheme, we include only a 3 _ in- 
elastic state of le O with f3 = 1.81 and excitation energy 
E* == 0.25 MeV. As will be shown in next section, this 
coupling is sufficient to reproduce the fusion cross sec- 
tion and measured barrier distribution. It is interesting 
to note that while both BPM and CCBPM predict a sat- 
uration value, CCDIFFUSION model does not show any 
saturation point. The < I > decreases with decreasing 
energy. This aspect can be verified experimentally if pos- 
sible by extending the measurements to lower energies. 

Like L(E) and S(E), another pragmatic way to rep- 
resent fusion data is through it's second derivative given 
by the relation, 



D(E) = 



1 d\a f E) 
*R\ 



dE 2 



(22) 



which is interpreted as the experimental barrier distri- 
bution smeared by a tunneling or transmission factor. 
Assuming T t (E) =T (e) where e = E -1(1 + l)h 2 /(^Rl), 
it can be shown that 



T (E) 



1 d(aE) 



nR 2 



dE 



(23) 



10" 
10 3 

10 2 
10 : 
10° 

CO 1° 2 

10 3 
1x10 4 
1x10 5 
10 6 
10"" 

40 

35 

30 

25 

A 20 

v 15 

10 

5 





t 58 Ni+ 58 Ni \ 64 Ni+ M Ni 



~i — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — r 




• DATA 

--CCBPM 

CCDIFFUSION 



_i , 1 1 i_ 



85 90 95 100 105 110 115 120 125 130 135 140 



1 I 1 I 1 I 1 


I 1 I 1 I 1 I 1 I 1 


16— 144— S 

0+ Smx 


f- CCBPM '_ 




f; CCDIFFUSION - 




BPM 




/ .' 58m- 58... 

/ .■ Ni+ Ni 


, i . r^T i , 


i . i , i . i . i , 



50 60 70 80 90 100 110 120 130 140 

E (MeV) 



FIG. 7: The S factor as a function of E cm (upper panel) for 
the same systems as shown in figure[6] For 64 Ni+ 64 Ni, the X- 
axis is shifted by +20 MeV for clarity. The lower panel shows 
< L > as a function of E cm for ie O + 144 Sm and 58 Ni + 58 Ni 
systems. The potential parameters are listed in table I. For 
58 Ni + 58 Ni, the X-axis is shifted by -10 MeV for clarity. 



Therefore, in case of a single barrier, the fusion cross 
section can be obtained by integrating the above equa- 
tion form to E and using Tq(E) either from Eq. (fT6|) or 
Eq. (fl"2|) . In the first case, 



a f (E,V b ) = ^log 



2tt 

^e XP i-(E- 



V h 



(24) 



where as in the second case, 



a f (E,V b )=nRl 



where 
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EV2n 



[XV^ (1 + erf X) + exp (-X)](25) 



X 



E-V b 



(26) 



and dvT = \/27 is assumed as constant. Note that while 
both Eq. ([24|) and Eq. ((25|) are the expressions for one di- 
mensional fusion cross section models, Eq. (|24p represents 
the celebrated Wong's formulae under parabolic approx- 
imation (BPM ) [25| where as Eq. ([25l) is identical to the 
expression used by Wilczynska and Wilczynski (See Eq. 
6 of Ref . ) to fit fusion cross section for several heavy 
ion systems. 
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FIG. 8: The D(E) versus E cm for 16 0+ 144 Sm and 58 Ni+ 58 
Ni systems. The data points are taken from Ref. The 
X-axis is shifted by -10 MeV for 58 Ni + 58 Ni system for 
clarity. 



The fusion cross section in Eq. (f21j) is written as the 
sum over the discrete barriers. Replacing this sum by 
an integral i.e. replacing the weight uo a by a distribution 
of barriers with probablity D(V) and using fusion cross 
section from Eq. (|24]) . it can be shown that 



where 



and 



1 d 2 (a f E) 
ttR 2 dE 2 



D(V)G(f)dV, 



expjf) 
[l + exp(f)Y 



f 



(27) 



(28) 



(29) 



If G(f), which is a sharply peaking function at E = V, 
can be replaced by a delta function S(E — V), Eq. ([27|) 
represents a true barrier distrbution as given by Eq. (|22[) . 
However, the above interpretation breaks down if we re- 
place T (E) by Eq. (fT2"|) . In case of diffusion, 



G(f) 



dE 



1 



T 



exp 



{E-Vf 
2 7 2 



(30) 



which is a Gaussian instead of a delta function and it will 
not be appropriate to interprete Eq. (|27|) as barrier distri- 
bution. Since all experimental measurements derived us- 
ing Eq. (|22[) are associated with true barrier distribution, 



it should be interpreted with some amount of caution. 
On the other hand, whatever be the interpretation, like 
L{E) and S(E), D(E) is another meaningful representa- 
tion of experimental data which is very sensitive to the 
coupling parameters (hence to the nuclear structure ef- 
fects) at and above the barrier region. Fig. [8]shows D(E) 
for 16 O + 144 Sm and 58 Ni + 58 Ni systems with the same 
coupling parameters as discussed before. The experimen- 
tal barrier distributions are shown only for 16 O + 144 Sm 
system. As can be seen, the predictions of both CCBPM 
and CCDIFFUSION are nearly identical and also both 
show the nuclear structure effects at and above barrier 
energies. Although er/ at deep sub-barrier energy is dif- 
ferent for both CCBPM and CCDIFFUSION , the D{E) 
is not sensitive to this varitioan. It may be mentioned 
here that unlike D(E), both L(E) and S(E) strongly de- 
pend on the variation of fusion cross sections at deep sub- 
barrier energeis. Therefore, any model, to remain valid 
from above barrier to deep sub-barrier energies, it should 
account for the above three factors simultaneously. 

Finally, we would like to mention here that although 
Eq. (|25|) is identical to the expression used in Ref. [IB] , the 
interpretations under which it has been derived are quite 
different. For example, the authors of Ref. (2(| start with 
a classical expression for fusion, 



^ = ^(1-1), 



which gives, 



1 d 2 (Ea) 
~dW 



nR 2 



= S(V b -E). 



(31) 



(32) 



Eq. ([231 can be obtained by replacing this delta function 
by a Gaussian distribution of barrier, 



P(V) 



1 



7V27T 



exp 



(E-v b y 

2 7 2 



(33) 



and integrating from to E. Using Eq. (|2l)]) and with 
Vb, Rb and 7 as parameters, the authors of Ref.[2(| have 
fitted fusion cross sections for several systems over a wide 
range of energies. Although, they are able to fit the data 
for several systems, the parameters turns out unphysical. 
For example, the Rb value for 58 iVi + 58 Ni is found out 
to be 6.0 fm where as the normal value is about 10.5 fm. 
Since the fusion cross section approaches its geometrical 
limit at higher energies (<7/ = ttR 2 ), the reduced radius 
parameter would result in a cross section which will be 
about 30% less. This unphysical behaviour is probably 
due to forcing an one diemnsional model to explain the 
fusion data over a wide range without using any channel 
coupling effects. 



IV. FUSION OF RADIO ACTIVE IONS 

In recent years, much efforts have been devoted to 
study the fusion reaction of light halo nuclei with stable 
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FIG. 9: The fusion cross section versus E cm for 6 He + 238 U 
and 6 He + 209 Bi systems and their corresponding S factors. 
The potential parameters are listed in table II. The experi- 
mental data points for 6 He -\- 209 Bi system are taken from 
Ref . . For 6 He + 238 U, the fission and fusion data points 
(squares and open circles) are taken from Ref. [14I] and the 
solid circles (normalized by a factor of 0.54 as suggested in 
Ref. [3) are taken from Ref.[T3|. The X-axis for 6 He + 209 Bi 
system is shifted by —10 MeV for clarity. 



heavy targets at energy around the Coulomb barrier. In 
light exotic nuclei such as 6 He, n Li and 11 Be, the neu- 
tron matter extends to a larger distance than the usual 
nuclear interaction scale. The fusion cross section for 
such system is expected to be enhanced with respect to 
the fusion with stable projectile due to the reduction in 
Coulomb barrier. On the other hand, the weak binding 
of the halo neutrons lead to a higher break up proba- 
bility of the projectile, which in turn, reduces the fusion 
cross section. Therefore, understanding the effect of neu- 
tron halo on fusion has been highly contoversial. In one 
of the earlier measurements in 6 He + 209 Bi reaction, a 
significance enhancement in fusion cross section with re- 
spect to 4 He + 209 Bi system was observed, espec ially 
at sub-barrier region [1JJ]. In the subsequent study 
an isolated 4 He group was observed with cross section 
far exceeding the fusion yield both above and below the 
barrier. This reaction mechanism was interpreted as the 
direct break up and transfer of neutrons to the unbound 
states in 211 BLa mechanism having close resembalnce to 
neutron flow [171 ] . The effect of the halo on the fusion- 
fission probabilty was also measured for 6 He + 238 JJ sys- 



tem [13j. The data lead to the conclusion that there is 
no hindrance to fusion at the barrier due to break up 
reaction and below barrier, there is a large enhancement 
of the fusion probability for the 6 He nucleus with re- 
spect to stable 4 He (see filled circles and squares in Fig. 
[9]). However, this measureemnt was confronted later on 
as no substantial fusion enhancement was observed [lj], 
rather, the large fusion yield was found due to two neu- 
tron tranfer induced fission of 240 U (The open circles in 
Fig. [5] which are obtained after correcting for transfer in- 
duced fission). It was also concluded that a simple BPM 
model may be adequate to explain the data. Recently, 
the yield of 210 Po isotope has been measured through 
the 2n evaporation channel of the fusion reaction e He 
with the 208 Pb target The yield demonstrates an 

extremely large enhancement of sub-barrier fusion cross 
section. Such a controversial situation in the experimen- 
tal data on sub-barrier reaction as it has been in the 
theoretical interpretations. 

Although we are still far from good understanding of 
the subject and much more experimental and theoretical 
studies are needed, in the following we examine the pos- 
sibility of explaining the fusion reaction involving halo 
nuclei based on the same diffusion model which has been 
used for heavy ion systems. Due to the extended neutron 
distribution, it is more probable that the system forms a 
di-nuclear complex at a large distance due to sequential 
neutron flow promoting neck formation. Thus, instead of 
tunneling, the di-nuclear configuration will diffuse over 
an adiabatic barrier in order to achieve fusion. Since the 
projectile is light, there should not be any coupling effect 
except the neck formation. The fusion process should 
follow an one dimensional diffusion model as given by 
Eq. (fT7j) with diffusion probability as given by Eq. (fT2|) . 
Figj9] shows (T/ as a function of E cm for 6 He + 238 U and 
6 He -f 209 Bi systems. The dotted curves show the pre- 
dictions based on the BPM model which has three pa- 
rameters, Vb, Rb, fiLj. The diffusion model uses the same 
Vb and Rb parameters, but the d is adjusted to get a rea- 
sonably good fit to the experimental data points. The 
lower panel shows the expected S factors. 



V. CONCLUSION 

In conclusion, we have proposed a new model based on 
the diffusion over an adibatic fusion trajectory to explain 
fusion cross section a/ and associated derived quantities 
like L(E), S(E) and D(E) for several systems both at 
above and below barrier energies. For the fusion of light 
projectile with heavy target for which neck fromation or 
neutron flow is not very meaningful (like fusion of He 
on stable neuclei), the system will fuse through one di- 
mensional tunneling of the sudden barrier and the fusion 
cross section can be described by Wong type calculation 
(a simple BPM model). However, if the projectile is a 
halo nuclei, but still lighter enough so that coupling ef- 
fects can be neglected (as coupling potential depends on 
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Z\Z-i), the fusion need to proceed through diffusion. Sim- diffusion as the mechanism of transmission over the adi- 
ilarly, for heavy projectile target combinations, coupled abatic Coulomb barrier, 
channel calculations are required to be carried out with 
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